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! ABSTRACT 

We present medium- resolution VLT/F0RS2 spectroscopy of six cataclysmic vari- 
ables discovered by the Sloan Digital Sky Survey. We determine orbital periods 
for SDSSJ023322. 61-1-005059.5 (96.08 ± 0.09 min), SDSS J091127.36-h084140.7 
qT ■ (295.74 ± 0.22min), SDSS J103533.02-h055158.3 (82.10 ± 0.09min), and 

^ ! SDSS J121607.034-052013.9 (most likely 98.82 ± O.lGmin, but the one-day aliases at 

• 92 min and 107 min are also possible) using radial velocities measured from their Ha 

' and H/3 emission lines. Three of the four orbital periods measured here are close to 

■ the observed 75-80 min minimum period for cataclysmic variables, indicating that the 

■ properties of the population of these objects discovered by the SDSS are substantially 

different to those of the cataclysmic variables found by other means. Additional 
photometry of SDSS J023322. 61-1-005059. 5 reveals a periodicity of approximately 

I ' 60 min which we interpret as the spin period of the white dwarf, suggesting that this 

O I system is an intermediate polar with a low accretion rate. SDSS J103533. 02-1-055158. 3 

' has a period right at the observed minimum value, a spectrum dominated by the cool 

^ ' white dwarf primary star and exhibits deep eclipses, so is an excellent candidate for 

an accurate determination of the parameters of the system. The spectroscopic orbit 
of SDSS J121607.03-f 052013.9 has a velocity amplitude of only 13.8 ± 1.6kms~\ 
k> \ implying that this system has an extreme mass ratio. From several physical con- 

. straints we find that this object must contain either a high-mass white dwarf or a 

' brown-dwarf-mass secondary component or both. 

Key words: stars: novae, cataclysmic variables - stars: binaries: close - stars: bina- 
ries: eclipsing - stars: white dwarfs - stars: individual: SDSS J023322.6H-005059.5 
- stars: individual: SDSS J091127.36+084140.7 - stars: individual: SDSS 
J093238.21-h010902.5 - stars: individual: SDSS J101037.05-H024914.9 - stars: indi- 
vidual: SDSS J103533.02+055158.3 - stars: individual: SDSS J121607.03-I-052013.9 - 
stars: individual: SDSS J133941. 11-^484727.5 



1 INTRODUCTION 

Cataclysmic variables (CVs) are interacting binary stars 
containing a white dwarf primary star and a low-mass sec- 
ondary star in a close orbit. In the majority of systems the 
secondary star is unevolved, fills its Roche lobe and loses 

* E-mail: j.k.taylor@warwick.ac.uk (JS), 
Boris . Gaensicke@warwick. ac.uk (BTG) , 
T.R.Marsli@warwick.ac.uk (TRM) 



mass to the primary star via an accretion disc. Comprehen- 
sive rev i ews o f the proper t ies of CVs have been given by 
IWarneJ lll995h and lHellieil feOOll) . 

The evolution of CVs depends primarily on angular mo- 
mentum loss (AML) from their orbit. CVs with long periods 
contain a secondary star which is only partially convective. 
It is thought that the presence of a magneti c field in this star 
causes AML by magnetic braking ( Vcrbun t fc Zwaanlll98ll : 
iRappaport. Joss. &i Webbinklll982il until the orbital period 
decreases to about three hours. At this point the secondary 
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star becomes fully convective and magnetic braking ceases. 
This change in the state of the star causes it to relax to its 
thermal equilibrium radius and so shrink inside its Roche 
lobe. Mass transfer ceases and weaker AML is now domi- 
nated by gravitational radiation JPaczviiskJll96'it) . Once the 
period reaches about two hours, the shrinking Roche lobe 
re-establishes contact with the secondary star, restarting 
the mass transfer. The period continues to decrease to an 
observed minimum of approximately 80 minutes, at which 
point the secondary star becomes a degenerate brown-dwarf 
like object. As the mass donor is now de generate, furthe r 
mass transfer causes the period to increase (|Pattersonll99jtl . 

Systems beyond this "period bounce" are theoreti- 
cally predicted to represent up to 70% of the steady- 
state population of CVs; population synthesis models 
also predict that the vast majority (99%) of C Vs have 
evolv e d to periods shorter than the 2-3 hour gap ijde Kooll 
Il992l : Ide Kool fc Ritteil Il993t iKolbl Il993l: IPolitanol Il99a 



)1 fc Ritt( 

313" 



iKolb fc Baraffelll999f) . The models also predict that the low 
rate of change of period near the minimum value causes an 
increase in the population of CVs there, which should man- 
ifest itself as a spike in the period distribution. Ifowever, 
whilst several candidates exist, not one CV has yet been con- 
firmed to have a secondary co mponent with a mass appro- 
priatc fo r a period bounce CV (IL ittlcfair. DhiU on. fc Martini 
[2003: Patterson. Thorstensen. fc Kem p 2005b). In addition, 
the observed period distribution of CVs does not have 
a spike at the minimum period and it does not contain 
far more short-period systems than those with periods 
beyond the well-known 2-3 hour gap (e.g. iDownes et alJ 
1200 ll: iRitter fc Kolbl l2003ll . Furthermore, the theoreti- 
cal minimum period is generally found to be around 
65 minutes, which is substantially shorter than the ob- 
served value. Many different AML mechanisms have been 
put forward as possible solutions to these problems 
fe.g.lAndronov. Pinsonneault. fc Sillsll2003l : pTaam fc SpruitI 
I2OOII: ISchenker et alJl200^ without complete success. The 
AML rates given by different prescriptions can differ by sev- 
eral orders of magnitude (^Schrciber fc Gansickc 200^. 

Due to these difflculties, detailed observational popu- 
lation studies are of great importance to provide the con- 
straints necessary to further develop the theory of compact 
binary evolution. The currently known population of CVs 
is afilicted by several strong biases which makes the charac- 
terisation of the true population - and a comparison with 
models - ex tremely difficult. The Sloan Digital Sky Sur- 
vey (SDSS^: lYork et alJ 1200(3) is a large project to obtain 
five-colour imaging and multi-object spectroscopy of a large 
area at high galactic latitudes in the Northern sky, with 
the main science goal being the measurement of redshifts 
for an extensive sample of galaxies and quasars. Objects 
for spectroscopic follow-up are chosen, on the basis of their 
photometric colours, to be blue and/or distant from the lo- 
cus of the stellar main sequence. Approximately 200 CVs, 
with a wide range of photometric colours, have been de- 
tcctcd in the SDSS spectroscopic database iSzkodv et alJ 
[2002, 2003, 2004, 2005, 2006), and many more are expected 
to be found from future SDSS observations. The limiting 
magnitude of this population is g ^ 20, much fainter than 
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previous surveys (e.g. ICreen. Schmidt, fc Liebertlll986l) . so 

the SDSS CV population is both more homogeneous and on 
average fainter than the currently known sample of CVs. It 
is therefore expected to contain more short-period systems, 
and may include a number of intrinsically faint CVs which 
have evolved beyond the period bounce. 

We are conducting a research program to meas ure the 
orbit al periods of SDSS CVs t.Gansicke,2005.; Gansicke et alJ 
I2OO6I : Dillon et al., 2006, in preparation) in order to assess 
whether the study of a less biased sample of CVs can 
reconcile observations with theories of the evolution of 
interacting binary stars. In this work we present phase- 
resolved spectroscopy of six faint SDSS CVs (Table 
and measure orbital periods for SDSS J023322.6-I-005059.5, 
SDSS J091127.36-h084140.7, SDSS J103533.02+055158.3 
and SDSSJ121607.03+052013.9. The spectra of 
SDSS J103533.02-I-055158.3 show that it is an eclipsing sys- 
tem. In addition we observed SDSS J101037.05-H024914.9, 
finding that it was several magnitudes fainter than in the 
SDSS observations, and SDSS J093238. 2-^010902. 5, but did 
not detect a coherent radial velocity motion. Our selection 
policy was to pick faint CVs which were discovered by 
the SDSS and were observable from the VLT. We did 
not attempt to carefully select a representative sample 
of objects because the goal of our research project is to 
(eventually) measure the orbital periods of all of the SDSS 
CVs with unknown periods. In this work we shall abbrevi- 
ate the names of the targets to SDSSJ0233, SDSSJ0911, 
SDSSJ0935, SDSS JlOlO, SDSSJ1035 and SDSSJ1216. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 VLT spectroscopy 

Spectroscopic observations were carried out in 2006 January 
using the FORS2 spectrograph on Unit Telescope 1 of the 
Very Large Telescope (VLT) at ESO Paranal, Chile. The 
1400V grism was used for the observations of SDSS J0233 
and the 1200R grism for observations of the other targets. 
Observations using the 1400V grism cover the wavelength 
interval 4630 to 5930 A, which includes H/3, at a recipro- 
cal dispersion of 0.64 Apx~^. Those using the 1200R grism 
were targeted at Ha and cover the wavelength interval 5870 
to 7370 A at 0.73 A px~^. A log of observations is given in 
Table 121 From measurements of the full widths at half maxi- 
mum (FWHMs) of arc-lamp and night-sky spectral emission 
lines, we estimate that the resolution of our observations is 

1.2 A at II/3 and 1.5 A at Ha. 

Data reduction was undertaken using optimal extrac- 
tion (.Home. .1986) as implemented in the PAMELA^ code 
jMarshlll98'9l) . which also makes use of the STARLINK^ pack- 
ages FIGARO and KAPPA. 

The wavelength calibration of the spectra was under- 
taken using the molly^ program and arc lamp exposures 
taken during daytime as part of the standard calibration 
routines for F0RS2. However, we have found that flexure 



^ PAMELA and MOLLY were written by TRM and can be found at 
http : //quetzel . esc .Warwick. ac .uk/phsaap/ software/ 



The Starlink Software Group homepage 
[http : //www. starlink.rl . ac.uk/. 
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Table 1. Apparent magnitudes of our targets in the SDSS ugriz passbands. r'gpcc is an apparent magnitude calculated by the SDSS 
by convolving their flux-calibrated spectra with the r passband function, and is obtained at a different epoch to the ugriz magnitudes 
measured from the CCD imaging observations. 



SDSS name 


Short name 


Reference 




u 


a 


r 


i 


z 


spec 


SDSS J023322. 61+005059.5 


SDSS J0233 


Szko^^^aL 


(2002) 


19.76 


19.86 


19.73 


19.93 


19.82 


19.92 


SDSS J091127.36+084140.7 


SDSS J0911 


Szko^^e^aL 


(2005) 


19.63 


19.70 


19.16 


18.66 


18.28 


19.51 


SDSS J093238. 21+010902.5 


SDSS J0932 


Szkoc^^^aL 


(2003) 


19.70 


20.26 


19.53 


19.10 


19.12 


19.71 


SDSS J101037.05+024914.9 


SDSS JlOlO 


Szkoc^^^aL 


' (2003'! 


20.39 


20.75 


20.35 


20.71 


20.46 


20.63 


SDSS J103533. 02+055158.3 


SDSS J1035 


Szkody et al. 


(2006) 


19.14 


18.79 


18.81 


18.97 


19.12 


18.99 


SDSS J121607.03+052013.9 


SDSS J1216 


Szkodv et al. 


(2004) 


20.32 


20.13 


19.98 


20.26 


20.27 


20.10 



Table 2. Log of the observations presented in this work. The acquisition magnitudes were measured from the VLT/FORS2 acquisition 
images and are discussed in Section l2.2l 



Target 


Date 


Start time 


End time 


Telescope and 


Optical 


Number of 


Exposure 


Acquisition 




(UT) 


(UT) 


(UT) 


instrument 


element 


observations 


time (s) 


magnitude 


SDSS J0233 


2004 08 25 


03 22 


05 50 


LT RATCAM 


SDSS g filter 


101 


60 




SDSS J0233 


2004 08 28 


03 41 


05 55 


LT RATCAM 


SDSS g filter 


92 


60 




SDSS J0233 


2005 01 02 


20 53 


23 55 


INT WFC 


SDSS g filter 


80 


60-120 




SDSS J0233 


2005 01 03 


20 04 


23 41 


INT WFC 


SDSS g filter 


113 


70 




SDSS J0233 


2005 01 04 


19 54 


22 54 


INT WFC 


SDSS g filter 


106 


50-80 




SDSS J0233 


2005 01 04 


20 38 


22 45 


WHT ISIS 


R600B R316R 


9 


900 




SDSS J0233 


2005 01 06 


19 59 


23 48 


WHT ISIS 


R600B R316R 


10 


900 




SDSS J0233 


2005 01 07 


00 08 


23 53 


WHT ISIS 


R600B R316R 


15 


900 




SDSS J0233 


2006 01 26 


00 36 


03 06 


VLT FORS2 


1400V grism 


18 


480 


19.6 


SDSS J0233 


2006 01 27 


00 29 


01 49 


VLT FORS2 


1400V grism 


10 


480 


19.9 


SDSS J0233 


2006 01 28 


02 18 


03 02 


VLT FORS2 


1400V grism 


10 


480 


19.9 


SDSS J0911 


2006 01 26 


04 22 


07 04 


VLT FORS2 


1200R grism 


19 


480 


18.4 


SDSS J0911 


2006 01 27 


03 40 


05 18 


VLT FORS2 


1200R grism 


12 


480 


18.4 


SDSS J0911 


2006 01 27 


06 54 


08 22 


VLT FORS2 


1200R grism 


11 


480 


18.4 


SDSS J0932 


2006 01 27 


02 09 


03 21 


VLT FORS2 


1200R grism 


9 


480 


19.9 


SDSS JlOlO 


2006 01 26 


03 26 


04 02 


VLT FORS2 


1200R grism 


5 


480 


22.3 


SDSS JlOlO 


2006 01 26 


07 17 


07 18 


VLT FORS2 


B filter 


1 


60 


22.0 


SDSS JlOlO 


2006 01 26 


07 19 


07 20 


VLT FORS2 


V filter 


1 


60 


21.9 


SDSS J1035 


2006 01 27 


05 40 


07 08 


VLT FORS2 


1200R grism 


22 


200 


18.8 


SDSS J1035 


2006 01 28 


03 22 


04 40 


VLT FORS2 


1200R grism 


20 


200 


18.7 


SDSS J1035 


2006 01 28 


08 39 


09 40 


VLT FORS2 


1200R grism 


16 


200 


18.7 


SDSS J1216 


2006 01 26 


07 32 


09 35 


VLT FORS2 


1200R grism 


14 


480 


19.5 


SDSS J1216 


2006 01 27 


07 22 


09 38 


VLT FORS2 


1200R grism 


16 


480 


19.7 


SDSS J1216 


2006 01 28 


04 56 


06 38 


VLT FORS2 


1200R grism 


11 


480 


19.8 



of the spectrograph can cause the position of the spectrum 
on the CCD to vary by up to 0.6 pixels (about 20kms^^) 
depending on elevation. We have removed this trend by ap- 
plying pixel shifts to the wavelength calibrations of the spec- 
tra to force the centroids of the night sky emissio n lines at 
5577.338 and 6300.304 A JOsterbrock et alJll996l) to these 
wavelengths. We have been able to check this procedure in 
some observations when other night sky emission lines are 
strong, and find that it is always accurate to within 0.1 pix- 
els. 



can be obtained. Exposure times were generally 20 s and the 
observations were unfiltered, which means that the Poisson 
noise is in all cases smaller than the systematic errors. We 
have extracted differential photometry from these images 
using the iraf* tool imexamine, taking special care to select 
comparison stars from the SDSS database with colours as 
close to our target CVs as possible in order to minimise 
colour effects for those acquisition images taken in white 
light. The resulting magnitudes (Table I^J approximate the 
SDSS r passband, with some uncertainty due to the fact 



2.2 VLT photometry 

The spectroscopic observing procedure of F0RS2 included 
obtaining target acquisition images from which photometry 



^ IRAF is distributed by the National Optical Astronomical Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc. under contract with the National 
Science Foundation. 
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that CVs have quite difFerent spectral energy distributions 
from the comparison stars. 

Images in the Johnson B and V passbands were also 
obtained for SDSS JlOlO after we found that spectroscopic 
observations of it had a much lower signal than expected. We 
used the colour transformations by R. Lupton on the SDSS 
website^ to convert the ugri magnitudes of the comparison 
stars to the B and V magnitudes needed for our differential 
photometry. 



2.3 WHT spectroscopy 

We obtained time-resolved spectroscopy of SDSS J0233 at 
the 4 m William Iferschel Telescope (WHT) on La Palma in 
2005 January, using the double-arm spectrograph ISIS with 
the R600B grating and 4kx2k pixel EEV detector in the 
blue arm, the R316R grating and 4.5kx2k pixel Marconi 
detector in the red arm, and a 1.2" slit width (Tabled . The 
blue arm spectra cover the wavelength range 3900-5400 A 
with a resolution of 1.8 A (corresponding to two binned pix- 
els) and the red arm spectra cover 6100-8900 A with a res- 
olution of 3.6 A. Regular arc lamp and flat-field exposures 
were obtained to adjust the wavelength calibration for the 
instrument flexure and to remove CCD fringing in the red 
spectra. The data were bias and flat-field corrected using 
FIGARO. Optimal extraction and wavelength and flux cali- 
bration were performed in PAMELA and molly. Due to the 
bad atmospheric conditions (variable seeing with an average 
of about 5 arcseconds) the quality of the individual spec- 
tra is relatively poor. The average spectrum of SDSS J0233 
from these observations shows strong double-peaked Balmer 
emission lines and relatively weak He I lines. 



2.4 INT and LT photometry 

Differential p-band photometry of SDSS J0233 was obtained 
on the Liverpool Telescope (LT) in August 2004 and the 
Isaac Newton Telescope (INT) in January 2005 (Table I^J. 
The INT was equipped with the Wide Field Camera (WFC), 
an array of four EEV 2kx4k pixel CCDs, and the LT with 
RATCAM, which contains an EEV 2kx2k pixel CCD. No 
binning or windowing was ado pted. The data were reduced 
using the pipeline described bv lCansicke et al] (1200411 . which 
performs bias and flat-field corrections within midas^ and 
uses the SEXTRACTOR package (Berlin & Arnouts 1996) to 
perform aperture photometry for all objects in the field of 
view. Instrumental g magnitudes were converted into ap- 
parent magnitudes using several comparison stars near the 
target along with their magnitudes given on the SDSS web- 
site. 



3 DATA ANALYSIS 

3.1 Radial velocity measurement 

We measured radial velocities from the hydrogen emission 
lines by cross-correlation with single and double Gaussian 
functions, as implemented in molly^. The single- Gaussian 
method gives radial velocities strongly affected by the peak 
of the emission line, so can highlight the velocity variation 
of the line emission of the bright spot on the edge of the 
accretion disc caused by the impact of material lost from 
the secondary star. 

T he double-Gaussian technique iSchneider fc Yound 
Il980l) can be used to measure the velocity variation of the 
wings of emission lines. The flux in the line wings comes from 
material close to the inner edge of the accretion disc, so is 
generally interpreted as representing the motion of the white 
dwarf. While common, th is interpretation is subject to sys - 
tematic uncertainties (see IShafteilll983 : lThorstensenll2000f) . 
The velocities are measured by cross-correlating the emis- 
sion line spectrum against a function containing two Gaus- 
sians of equal FWHM and area but with one having positive 
height and one negative height. By choosing a small FWHM 
(similar to the resolution of the spectra), and a separation 
between the two Gaussians of slightly less than the FWHM 
of the emission line, the point at which the cross-correlation 
function is equal to zero represents the velocity of the line 
wings. In all cases we varied the FWHM and separation of 
the Gaussians to confirm the robustness of our results. 

An independent orbital period determination can be 
carried out using the ratio of the fiuxes in the redder and 
bluer parts of the line. Our procedure for this is to deter- 
mine the wavelength where the line fiuxes bluewards and 
redwards are equal, expressed as a velocity relative to the 
rest wavelength of the line. When reliable results were ob- 
tained from this method, they agreed in all cases with those 
we found from the radial velocities. 

Averaged emission line profiles are given in Fig. Q and 
the SDSS spectra of the CVs for which we have measured 
orbital periods are shown in Fig.|5|for reference. 

3.2 Orbital period measurement 

The measured radial velocities for each CV were 
searched for peri o ds us ing periodograms computed us- 
ing the IScargl el il982^). analvsis of variance (AoV; 
LSchwarzcnbcrg-Czcrnv '1989") and orthogonal polynomial 
(ORT; .Schwarzenberg-C zernv 199 6) methods implemented 
within the TSA* context in MIDAS. Periodograms for our four 
main targets are shown in Fig.|3 

We fitted single-lined circular spectroscopic orbits 
(equivalent to fitting sine curves) to the radial velocities us- 
ing the SBOP^ program, which we have previously found to 

! ive reliable error estim ates for the optimised parameters 
Southworth et aljlioosi) . 



^ Equations are given relating ugri to B and V, and are 
attributed to Lupton (2005) without further information about 
this reference. For further information see 

[http : //www . sdss . org/dr4/ algorithms / sds sUBVRI Transf orm . html | 
" ,http : //www. eso . org/projects/esomidas/. 



The reduced spectra and radial velocity observa- 
tions presented in this work will be made available 
at the CDS I http://cdsweb.u-strasbg.fr/ 1 and at 
http : //www. astro . keele . ac .uk/^^jkt/ 

^ http: //www, eso , org/projects/esomidas/ doc/user /98NQV/volb/node 220 . html 

^ Spectroscopic Binary Orbit Program, written by P. B. Etzel, 
[http : //mintaka. sdsu. edu/f aculty/etzel/ 
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SDSS J0233 SDSSJ0911 SDSS J0932 SDSSJ1010 SDSSJ1035 SDSSJ1216 
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Figure 1. The averaged emission line profiles of tlie six CVs for wfiich we obtained spectra (H/3 for SDSS J0233 and Ha for the others). 
The abscissa of each plot is in counts where the continuum level has been normalised to unity. 
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Due to the limited observing time available, we were not 
able to obtain data over a sufficient time or range of hour 
angles to avoid the periodograms containing strong power at 
the one-day aliases of the orbital periods of the CVs. In sev- 
eral cases this has made it difficult to be certain which of the 
possible peaks represents the true orbital period. When this 
problem has been encountered, we have fitted spectroscopic 
orbits to the radial velocities for each of the possible periods 
(Table In general, we expect the orbit with the correct 
period to have both the lowest residuals of the fit, as the sys- 
tematic deviations from the sine curve are minimised, and 
the highest velocity amplitude, as any mismatch in orbital 
phase would act to decrease the measured amplitude. 

As a final check, we have also investigated the reliabil- 
ity of our period determinations using b ootstrapping simu- 
lations, following the method outlined bvlPress et alJ lll992l . 
chapter 15.6). For each simulation we identified the point 
of highest power, which always fell close to one of the alias 
periods. The fraction of 'successes' for each alias can then be 
interpreted as the probability that it represent the correct 
period. We would like to emphasise that this fraction gener- 
ally represents an underestimate of the probability that the 
alias is the correct one. This is because it does not take into 
account the interactive rejection of aliases, normally because 
they result in clearly deformed radial velocity curves, which 
is possible when analysing the periodograms of the actual 
observational data for each system. 

For our four main targets we show the spectroscopic 
orbits in Fig.|l] phased according to our preferred ephemeris 
for each system, and give the orbital parameters in Table |1| 
The corresponding phase-binned spectra are shown in Fig.|S] 
(stacked) and Fig. |S| (trailed) . 



5000 6000 7000 8000 9000 

Wovelength (A) 



Figure 2. The SDSS spectra of our four main targets, all plot- 
ted using the same limits in wavelength and flux. The flux levels 
have been smoothed with a 10-pixel Savitsky-Golay filter using 
the SAVGOL procedure in IDL (see http://www.rsinc.com/ I. The 
units of the abscissa are 10~^^ Wm~-^ A~^, which corresponds to 
' erg s~ 



4 RESULTS FOR EACH SYSTEM 
4.1 SDSS J023322. 6+005059. 5 

Th e discovery spect rum of SDSS J0233 was presented 
by ISzkodv et al.l i2002i) and shows slightly double-peaked 
Balmer emission lines along w it h wea k He I emission (Fig.|5]l. 
IWoudt. Warner, fc Pretoriuj (|20^ obtained high-speed 
photometry of the system, finding both flaring and flickering 
activity but no clear indications of periodic behaviour. 

The higher quality of the WHT spectroscopy compared 
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Table 4. Circular spectroscopic orbits found using SBOP. The reference times (corresponding to zero phase) refer to inferior conjuction 
(the midpoint of primary eclipse in eclipsing systems). 



Target 


Orbital period 


Reference time 


Velocity amplitude 


Systemic velocity 


^rms 




(day) 


(HJD) 


( kms"-') 


(kms-l) 


( km S~^) 


SDSS J0233 


0.066725 ±0.000064 


53762.59652 ±0.00080 


44.8 ±3.4 


95.6±2.3 


12.7 


SDSS J0911 


0.205374 ±0.000152 


53762.61735 ±0.00062 


106.4 ±2.3 


53.3±1.6 


9.2 


SDSS J1035 


0.056983 ±0.000061 


53761.75152 ±0.00176 


60.0±3.7 


24.7±2.7 


18.2 


SDSS J1216 


0.068628 ±0.000112 


53762.81177 ±0.00123 


13.8 ±1.6 


7.3±1.1 


7.0 



SDSS J0233 H/3 SDSS J0911 Ha 




-2000 -1000 1000 2000 -1000 -500 500 1000 

Velocity (km s"') Velocity (km s"') 



Figure 5. Phase- binned stacked spectra of the target 



Table 3. Some spectroscopic orbital parameters for the possible 
periods for SDSS J0233, SDSS J0911 and SDSS J1216. ffrms is the 
root mean square of the residuals of the observations around the 
best fit. The adopted periods are shown using bold type. 



Target 


Period 


Velocity amplitude 






(day) 


( km s^-"-) 


( km s"-"-) 


SDSS J0233 


0.062692 


42.2 ±4.0 


15.2 


SDSS J0233 


0.066724 


44.8 ±3.4 


12.7 


SDSS J0233 


0.071446 


44.3 ±3.7 


13.7 


SDSS J0911 


0.256545 


124.4 ±4.0 


12.8 


SDSS J0911 


0.205374 


106.4±2.3 


9.2 


SDSS J0911 


0.171225 


98.3 ±4.7 


20.5 


SDSS J1216 


0.059963 


12.7±1.9 


8.1 


SDSS J1216 


0.064001 


13.5 ±1.7 


7.4 


SDSS J1216 


0.068628 


13.8±1.6 


7.0 


SDSS J1216 


0.073986 


13.7±1.7 


7.3 



to the SDSS data clearly shows that the Balmer lines are 
double-peaked, and suggests that SDSS J0233 is a disc- 
accreting system. We were not able to extract a sensible 
radial velocity variation from the strong Ha line, but a 
single-Gaussian approach to H/3 indicated an orbital period 
of ~ 90 min. We therefore decided to obtain the VLT spec- 
troscopy of SDSS J0233 centred on H/3. 

We have measured radial velocities from the H/3 emis- 
sion line in our VLT spectra using both single and double 
Gaussians. Radial velocities from the double Gaussian tech- 
nique show a weak periodicity at 131.4 and 120.6 minutes, 
but those from single Gaussian measurements show a clear 



SDSS J1035 Ha SDSS J1216 Ha 




-2000 -1000 1000 2000 -2000 -1000 1000 2000 

Velocity (km s"') Velocity (km s"') 



1. The top spectrum in each plot is the averaged spectrum. 

periodicity in the region of 96 min (Fig. I^J which results 
in a sinusoidal velocity variation with orbital phase (Fig.^. 
Trailed greyscale plots of the VLT spectra, phase-binned us- 
ing this period, show that a weak S-wave is present. S-waves 
are commonly observed in CVs and are attributed to the 
velocity variation of a hot spot where the accretion stream 
from the secondary star hits the outer edge of the accretion 
disc. Because of the larger amplitude of the S-wave, com- 
pared to the wings of the emission lines, it is often a better 
indicator of orbital variation than the emission line wings. 

The Scargle periodogram for SDSS J0233 (Fig. 01 and 
the CRT periodogram show the strongest power at 96 min, 
but the two adjacent one-day aliases show only slightly lower 
power. In the AoV periodogram the 90 min alias shows a 
much weaker peak but the 102 min peak again has only 
slightly less power than the 96 min highest peak. We have fit- 
ted spectroscopic orbits using all three possible periods and 
find that the 96 min period gives the orbit with the lowest 
residuals and the highest velocity ampitude. Bootstrapping 
simulations using the Scargle periodogram method find that 
this alias has the highest power 75% of the time (but remem- 
ber that this is pessimistic - see Sec. 13.21 . We are therefore 
confident that this is the actual orbital period of the binary 
system. The corresponding best-fitting orbit has a period of 
96.08 ±0.09 min. 

4.1.1 Light curve analysis for SDSS ,10233 

The light curves of SDSS J0233 display a mixture of undu- 
lations on timescales of hours, superimposed by short-term 
flares lasting about 60 min (Fig.|7J. These short-term bright- 
enings were also detected in the photometry reported by 
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Figure 6. Greyscale plots of the phase-binned trailed spectra of the target CVs. 



IWoudt et alj ||2004|) . Similar short events are seen on oc- 
casion as brightness depressions, e.g. on 2005 January 4 
(Fig. |7|l. A time series analysis of all data combined does 
not reveal any coherent period. However, different subsets 
of data show periodic variability. More specifically, a Scar- 
gle periodogram of the two nights of LT data obtained in 
2004 August (Table |^ shows a strong signal near 160 min 
(9 cycled"^); the data folded over this period show a com- 
plex but clearly repeated structure (Fig. |51 lower panels). 
On 2005 January 3 a strong ~60 min modulation is seen in 
the data and picked up by a Scargle periodogram (Fig. |H1 
top). There is no indication of this period in the spectra. 
Given the spectroscopically determined period of 96 min, we 
tentatively interpret the photometric signals as the white 
dwarf spin period of ~ 60 min and the beat between orbital 
and spin period of ~ 160 min. From this we suggest that 
SDSS J0233 is an intermediate polar with a low accretion 
rate. This is in agreement with the evid ence for strong X- 
ray variability reported bv lSzkodv et al.l l(2002'l . 



4.2 SDSS J091127.36+084140.7 

ISzkodv et al.l i2005ft presented the discovery of this CV on 
the basis of an SDSS spectrum which shows strong hydrogen 
Balmer emission and some Hei line emission superimposed 
on a slightly red continuum (Fig. |5J. 

We obtained a total of 42 Ha spectra from three nights 
of VLT observations. Radial velocities were measured using 
the double Gaussian technique with Gaussian FWHMs of 
300kms~^ and separation of 1500 kms~^. All three types of 
periodogram show strong power at periods close to 295 min. 
The one-day alias at 245 min has the strongest power in the 
Scargle periodogram, but phasing the radial velocities with 
this period gives a discontinuous radial velocity curve which 
suggests that the Scargle method does not perform as well 
when the phase coverage is incomplete. The AoV and ORT 
alternatives show only weak power at this period but power 
at the other one-day alias, 369 min, which is nearly as strong 
as for the 295 min period. Inspection of the phased data 
shows that the 295 min period gives nice sinusoidal radial 
velocity variations whereas the 369 min alias again result in 
a clearly deformed radial velocity curve. We are confident 
that an orbital period of close to five hours is correct, and 



fitting the radial velocities with SHOP gives a spectroscopic 
orbit with a period of 295.74 ± 0.22 min. 

At the suggestion of the referee, we have compared the 
red part of the SDSS spectrum of SDSS J0911 to spectral 
templates of M-dwarf stars in order to determine the spec- 
tral type of the secondary component. We find M2±0.5 from 
the relative strengths of the visible molecular absorption 
bands. Using the calibration of th e surface brightness in the 
7165 A and 7500 A TiO bands bv iBetiermann fc Weichholdl 
(^199?), and assuming a Roche-lobe-filling mass donor with 
radius R2 = (3.7±0.2) x 10* m, we estimate the distance to 
SDSS J0911 to be d = 2000 ± 300 pc. This distance and the 
galactic latitude of the system, b = -1-35°, imply that it is 
located ~ 1.1 kpc above the galactic plane. 



4.3 SDSS J093238. 2+010902. 5 

SDSS J0932 was discovered to be a CV bv ISzkodv et alJ 
l)2003h from its SDSS spectrum, which shows a faint flat 
continuum with reasonably strong hydrogen emission lines 
and slight Hei emission. 

We obtained a total of nine VLT spectra of SDSS J0932 
over an 80 min time period. These show Ha emission with a 
FWHM of about 900 km s"^ (Fig.IB, which is smaller than 
in the four systems for which we have measured radial ve- 
locity variations. The He i 6678 A line is also clearly in emis- 
sion, with an equivalent width of about one sixth of Ha. Our 
preliminary reductions at the telescope gave no clear indica- 
tion of radial velocity changes due to orbital motion, so no 
further observations were obtained. Fig. 1^1 shows radial ve- 
locities measured using the double Gaussian technique with 
a FWHM of 300 kms"^ and separation of 1000 kms"^; sim- 
ilar results are obtained for a range of values for Gaussian 
FWHM and separation. 

It is certainly possible that this is a system with an 
orbital period much longer than the length of our dataset, 
which was observed over an interval in its orbit when its 
radial velocity was changing by very little. If this is the case 
then our survey will be slightly biassed against longer-period 
CVs, so we intend to revisit SDSS J0932 at a later date to 
further investigate this system. 
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Figure 8. Time series analysis of the LT and INT photometry of SDSS J0233. Upper panels: The Scargle periodogram of the INT 
data of 2005 January 3. Left: a strong peak is present at 24 cycled"^ (OOmin). Right: the data folded over that frequency (small points) 
and binned into 15 points. Lower panels: the Scargle periodogram of the two combined nights of LT data from 2004 August. Left: the 
strongest signal is found at 8.7cycled~^ (166 min), with a number of aliases spaced by ^cycled"^. Right: light curves of the original 
data folded with 8.7cycled~^ (small points) and 25-point phase-binned data (large points). 



4.4 SDSS J101037.05+024914.9 

SDSS JlOlO was identified as a CV bv lSzkodv eTsI\ H 
from its SDSS spectrum, which shows a faint flat continuum 
with exceptionally strong hydrogen emission lines and faint 
He I line emission. 

The apparent magnitude of SDSS JlOlO in both its 
SDSS imaging and spectroscopic observations was around 
20.5, which is well within reach of VLT/FORS2 spectroscopy 
using our observational setup. However, the five spectra we 
obtained contained about ten times fewer counts than ex- 
pected, indicating that the object was much fainter than 
during the SDSS observations. The average spectrum is 
shown in Fig Q and displays weak Hq emission. 

Once photometry on our acquisition image confirmed 
that SDSSJIOIO was dimmer than expected, we obtained 
images in the B and V passbands. The magnitudes of 
SDSSJIOIO in these images is close to 22.0, confirming that 
SDSSJIOIO was in a lower state at the time of our observa- 
tions than those of the SDSS photometry and spectroscopy. 

We have measured the equivalent width of the Ha line 
to be approximately 600 A in the SDSS spectrum. This is ex- 
ceptionally large even for a short-period dwarf nova and sug- 
gests that the continuum is very weak, perhaps due to a cold 
white dwarf and a very late-type secondary star. The equiva- 
lent width of Ha in the average spectrum from our VLT data 
is only 43 ± 5 A, so the emission has weakened substantially 
in the interval between the SDSS and our VLT observa- 
tions. Further observations of SDSS JlOlO are desirable, as 
this object may be an example of the low-luminosity evolved 
CVs whose existence is predicted by theoretical models of 
CV evolution but which have not yet been observationally 
confirmed. 



4.5 SDSS J103533. 02+055158. 3 

SDSS J103 5 is a particularly interesting system; its SDSS 
spectrum dSzkodv et alJ l200(jl shows relatively weak but 
double-peaked Balmer emission fines superimposed on a 



very blue continuum with unmistakable wide Balmer ab- 
sorption lines coming from the white dwarf in the system 

(Fig.m. 

The brightness of SDSS J1035 allowed us to reduce our 
VLT/FORS2 exposure times to 200 s, resulting in a total 
of 58 spectra. These were measured for radial velocity, with 
best results obtained using double Gaussians separated by 
2000 km s"^ and with FWHMs of 300 km s"^ The disap- 
pearance of the blue or red emission line components in cer- 
tain observations (Figs. |K| and |SJ indicates that SDSS J1035 
is an eclipsing system. The changing line shapes through 
eclipse have strongly affected a few of the radial velocities, 
so we searched for periods using the AoV method. We find 
a period of 82.06 ± 0.09 min (Fig. from fitting the radial 
velocities with SHOP. There is an additional possible period 
at 65 min, where the two observed eclipses occur at the same 
phase, but this gives a clearly poorer radial velocity curve 
and so is very unlikely to refiect the actual period of the 
system. The 82 min period is confirmed by subsequent pho- 
tometry (Littlefair et al., submitted). For this period, the 
eclipses occur at phase 0.85 in the emission-line spectro- 
scopic orbit, which is quite distant from the expected posi- 
tion at phase zero for a Keplerian orbit. This efi'ect has been 
noted in quite a few CVs and indicates that radial velocities 
obtained from emission line wings are n ot always a good in - 
dicator of the motion of the white dwarf jThorstensenl2000ri . 



4-5.1 A model for the optical spectrum of SDSS J1035 



Followin g the approach d e scribe d in 
[iggi), iRodrfguez-Gil et al.1 ll2005l) and 



Gansicke et al 



Gansicke et al 



||2006D . we have fitted a three-component model to 
the SDSS spectrum of SDSSJ1035. The three compo- 
nents are a synthetic pure hydrogen white dwarf spec- 
tru m calculated with the tlusty and synspec codes 
( Hu benv fc Lan^ |l995l), an isothermal and isobari c hy- 
drogen disc jG^isickeT^Beuermann. fc Tho mas' "1997') and 
an observed late-type spectral template (M0.5-M9 from 
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Figure 3. Periodograms of the radial velocities measured for 
SDSSJ0233, SDSSJ0911, SDSSJ1035 and SDSSJ1216. That 
for S DSS J1035 is an OR T periodogram ISchwarzenberg-Czernvl 
Il996h and the others are IScargld ll982l) periodograms. The fre- 
quencies corresponding to the best-fitting orbital periods from 
SBOP are indicated with downward-pointing arrows. 



[Beuermann et al] Il998l and L0-L8 from iKirkpatrick et all 
Il999r) . We fixed the surface gravity of the white dwarf model 
spectra to \ogg = 8.0, corresponding to Mwd ~ 0.6 Mq, 
and a radius of iiwD ~ 8.7 x 10^ cm (assuming a carbon - 
core mass-radius relation from iHamada fc SalpeteJ Il96ll) . 
The radius of the secondary star is constrained by the size 
of the secondary Roche lobe, which primarily depends on 
the orbital period, and we fixed R2 ~ 8.5 x 10^ cm. Free pa- 
rameters of the three-component model are the white dwarf 
temperature, Two, the spectral type of the secondary star, 
SpT(2), the temperature Tdisc and the column density Sdisc 
of the hydrogen slab, and the distance, d, to the system used 
to scale all three components. 

A plausible model of the SDSS spectrum of SDSS J1035 
is achieved for the following set of parameters (Fig. llOU : 
TwD = 11 700 K, SpT(2) =L0, Tdi^c = 6500 K, Edisc = 5.3 x 



SDSS J0233 



150 - I 




0.0 0.2 0.4 0.6 0.8 1.0 

Orbitol phase 



Figure 4. Spectroscopic orbits for SDSSJ0233, SDSSJ0911, 
SDSSJ1035 and SDSS J1216. 



10~^gcm~^, and d — 280 pc. The radius of the accretion 
disc implied by the flux scaling factor, at a distance of 280 pc 
and assuming an inclination of 85° , is 7?disc « 2 x 10^" cm, 
well within the Roche lobe of the primary. The spectral type 
of the secondary star is constrained by the non-detection of 
TiO bands in the red part of the spectrum of SDSS J 1035 
and represents only an "early" limit - the true spectral type 
could well be later than LO. The white dwarf contribution to 
the optical spectrum of SDSS J1035 is approximately 85%. 

The parameters of the white dwarf are subject to un- 
certainty in its mass, which cannot be constrained from 
the optical spectrum alone, and to the contamination of 
the Balmer absorption lines by emission from the accretion 
disc. We estimate the uncertainties in the temperature of 
the white dwarf and distance to be approximately ±1500 K 
and ±60 pc, respectively. Adopting a larger surface gravity 
of log (7 = 9.0 for the white dwarf (corresponding to a mass 
of 1.2 M0) causes the best-fit parameters for the spectrum 
of SDSSJ1035 to change to Two = 14 000 K, d = 140 pc. 
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Figure 7. Sample light curves of SDSS J0233. 
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Figure 9. Measured Ha radial velocities of SDSS J0932. 



and an upper "early" limit on the secondary spectral type 
of L2. 

Comparing the properti es of SDSS J1035 to those of 
SDSS J133941. 11+484727.5 JCansicke et alJ l2006ll . it ap- 
pears that the systems are extremely similar, with orbital 
periods close to the observed 80 min period minimum and 
relatively cool white dwarfs which nevertheless dominate the 
optical spectrum of the systems. The implications are that 
both systems have very low mass transfer rates, and hence 
low-luminosity accretion discs and bright spots. They are 
good candiates for brown-dwarf-mass donor stars, as even 
the presence of a late-type M dwarf is difficult to reconcile 
with the non-detection of TiO absorption bands in the red 
part of the spectrum. Both systems thus have in all aspects 
the properties predicted by CV evolution theory for old CVs 
that have evolved to, or past, the minimum period. 




4000 5000 6000 7000 8000 

Wavelength [A] 



Figure 10. Three-component model of the SDSS spectrum of 
SDSS 1035. The three components are a white dwarf with T-yvD = 
11 TOOK, i?wD = 8.7 X 10* cm; and isothermal and isobaric hy- 
drogen slab with T^igc = 6500 K and T,^^sc = 5.3 X 10"-^ gcm~-^, 
and an LO secondary star, all scaled to a distance of d = 280 pc. 



4-5.2 Doppler tomography of the accretion disc of 
SDSS J 1035 

The prominent S-wave visi ble in SDSS J 1035 (Fig.EJ is com- 
parable to that of WZ Sge JSpruit fc Ru tton 1998). There is 
dynamical information in such S-waves which we attempted 
to extr act by application of the method of Doppler tomog- 
raphy jMarsh fc Hornelll988l) . In this procedure the spectra 
are modelled as the sum of a series of sinusoids which have 
radial velocities as a function of orbital phase (j> obeying the 
relation 

V{<j)) — —Vx cos 2-K(\> + W sin 27r0 

A Doppler map is a plot of the intensity of each sinusoid as 
a function of Vx a nd W . We a pplied the maximum entropy- 
based method of iMarsh fc H ornc (1988). A complication 
in the case of SDSS J1035 is the deep central absorption 
clearly visible in Fig. |S| This is lower than any emission 
pattern over the disc can produce, and is presumably the 
result of the underlying white dwarf photospheric absorp- 
tion line with possibly some contribution from absorption 
through the disc along the line of sight to the white dwarf. 
This causes a slight difficulty when applying the maximum 
entropy method, which assumes positive intensities at all 
points. Here we sidestep this problem by first adding Gaus- 
sians of FWHM eOOkms"^ to all the profiles of suflicient 
strength to allow the corresponding Doppler map to be pos- 
itive. We then computed the map (a two-dimensional Gaus- 
sian) corresponding to the Gaussians and subtracted this 
from the maps. This procedure does not affect the form of 
the map because we use an entropy that is only sensitive to 
structure on small scales (less than 250kms~^). 

The maps were computed avoiding phases from —0.08 
to +0.08 to avoid spectra affected by the eclipse, as eclipses 
are not accounted for in standard Doppler tomography. We 
used the ephemeris 

To (HJD) = 2453799.4807237 + 0.0570076 x E 
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Figure 11. The Doppler image in Ha of SDSS J1035. The Roche 
lobe of the mass donor and tracks of the gas stream for the param- 
eters indicated are shown. The lower track is the direct velocity 
of the stream while the upper represents the Keplerian velocity of 
the disc along the path of the stream. The left-hand panel shows 
the Doppler image whilst the right-hand panel shows the Doppler 
image after the subtraction of an azimuthal average to emphasise 
the asymmetries. 
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Figure 12. The Ha profiles of SDSS J1035 are shown on the left, 
with integer cycle gaps removed. The middle panel shows the fit 
based upon the map of Fig. 1111 while the right-hand panel shows 
the difference between the two. The greyscales of left-hand and 
middle panels run from to 5 while the right-hand panel is scaled 
from —2.5 to +2.5, in units of the continuum level. 



derived from independent photometric observations (Lit- 
tlefair et al., submitted). The map is shown in the left- 
hand panel of Fig. 1111 the right-hand panel shows the re- 
sult of subtraction of an azimuthally averaged version of 
the map to bring out the asymmetries. Fig. 1121 shows the 
spectra, fit and residuals. We plot a representative path of 
the gas stream that goes through the bright spot. Equally 
good fits can be obtained for many values of the mass ra- 
tio, q, and the sum of the velocity amplitudes of the two 
stars, (kms~^). Examples are (g,-ft'wD+J^2) = 

(0.02,550), (0.03,590), (0.05,620) and (0.07,650), although 
the last pair imply a white dwarf mass in excess of the Chan- 
drasekhar limit. It is not, however, certain that the bright- 
spot emission tracks the velocity of the gas stream; indeed 



the map here is reminiscent of the Doppler map of U Gem 
for which the emis sion is thought to lie betwe en the two 
tracks plotted here lIMarsh et alJlll99Q : .Unda-Sanzana et alJ 
2006). This would reduce the value of -ffwD+^f2 necessary to 
match the bright spot and hence relax the constraint upon 
the mass ratio implied by the Chandrasekhar limit. One fi- 
nal noteworthy point is the velocity of the ring from the 
outer disc visible in the map which peaks at a radius of 
about 800 km s"'^. This is one of the highest in our expe- 
rience; it suggests a combination of a fairly massive white 
dwarf and/or a small disc. 

4.6 SDSS J121607.03+052013.9 

SDSS J1216 was discovered to be a CV by ISzkodv et alJ 
(2004) from its SDSS spectrum, which shows a faint flat 
continuum with double-peaked hydrogen emission lines ex- 
hibiting a strong Balmer decrement. There is also faint He I 
emission which is strongest in the 5876 A line (Fig. |5J . 

We measured radial velocities of the Ha line in our 
41 VLT spectra using both double Gaussians (separation 
2000 km s"^ and FWHM 300 km s"^ provided the best re- 
sults) and a single Gaussian function (FWHM 2000 kms~^). 
The Scargle periodogram of the resulting radial velocities 
contains a strong signal at two aliases, with almost equal 
power, at 92min and 99min. The two adjacent aliases of 
86min and 107 min also have almost as strong power. The 
AoV periodogram clearly prefers the 99 min period and has 
only weak power at 106.5 min. The ORT periodogram, how- 
ever, prefers 92 min, with the second-highest power being 
at 99 min. On the basis of the periodograms alone, there- 
fore, we have identifled the 99 min period as most the likely 
to represent the orbital motion, but cannot fully rule out 
the 92 min period. Bootstrapping simulations confirm that 
99 min is the most likely period but that there is a signifi- 
cant probability that 92 min (30% for the Scargle and ORT 
periodograms, 10% for AoV) or 107 min (about 20% for the 
three periodigram types) is the correct one. 

We have fitted spectroscopic orbits to the radial ve- 
locities using the four possible periods, finding that the 
99 min alternative gives both the smallest residuals and 
the largest velocity amplitude, with the 107 min option 
marginally second-best in both criteria. As our methods 
generally prefer the 99 min period, with either the longer 
or shorter one-day aliases of this being the second-best op- 
tion, we expect that this period is the correct one. We have 
fitted the radial velocities with SBOP to find a period of 
98.82 ±0.16 min, resulting in a radial velocity curve which is 
sinusoidal. This period is also reproduced in studies of the 
ratio of the fiuxes in the blue and red parts of the emis- 
sion lines (see Sec. 13.111 . However, given that we have not 
been able to fully dismiss the alternative periods either side 
of 99 min, and that this object has other particularly inter- 
esting properties, a futher study of SDSS J1216 would be 
desirable. 

^.6.1 Nature of the donor star in SDSS J1216 

The spectroscopic orbit we have measured for SDSSJ1216 
(Table |3J has a very small velocity amplitude of 13.8 ± 
1.6kms~^, despite the presence of double-peaked line emis- 
sion suggesting a high inclination, indicating that the system 
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Figure 13. Diagrams of white dwarf mass versus orbital inclination (left) and secondary star mass (right) for SDSS J1216 (see Sec- 
tion l4.6.Tl for further explanation). The dark shaded regions indicate areas where inclination must be 75° or more (which is not expected). 
The light shaded regions show where an accretion disc on a Keplerian orbit at the observed velocity of the edges of the Ha emission 
line would be inside the surface of the white dwarf, which is not possible. The dashed lines indicate the edge of this region. The solid 
lines indicate the locus of points satisfying the mass function for orbital inclinations of 75°, 60°, 45° and 30° (labelled). The canonical 
minimum mass for core hydrogen burning is shown using a dotted line. 



has an extreme mass ratio. We can apply several constraints 
on the possible properties of the component stars in the sys- 
tem (see Fig. I13II :- 

• We have assumed in this discussion that the velocity 
amplitude found from the line wings is representative of the 
actual orbital motion of the white dwarf. This assumption 
is likely to be only approximately true, but the emission 
hne motion is generally accepted to overestimate the actual 
movement of the white dwarf, so this is a conservative as- 
sumption for the discussion below. It would be useful to 
measure the WD absorption features directly, which may be 
possible in the future from ultraviolet spectroscopy of the 
system. 

• A relation between masses and radii of white dwarfs was 
obtained (P. Bergeron, private communication) assuming 
Teff = 15 000 K, supplemented with the zero-temperature 
predictions of iHamada fc Salpeteil il96 j) from 1.2 M0 to 
1.4 M0. Our use of two different sources and our choice of 
Tcft makes a negligible difference to the results found below. 

• The full width at zero intensity of emission of the Ha 
hne gives a lower limit on the mass of the white dwarf. This 
quantity is notoriously t ricky to mea sure for wide emis- 
sion lines of faint stars JWarnejll995l) so we constructed 
a summed Hq spectrum in which the (very small) orbital 
motion of the line wings was removed. The full width was 
measured to be 100 ± 10 A, corresponding to a half- width 
zero- intensity velocity of 2290 ± 230kms~^. Assuming the 
mass-radius relation for white dwarfs we adopted above, this 
velocity gives the orbital inclination as a function of mass 
for which the innermost part of the disc is on a Keplerian 



orbit at the surface of the white dwarf. This provides an 
upper limit on the white dwarf radius which translates into 
a lower limit on its mass. 

• The spectra show no sign of eclipses (which can be seen 
by comparing the trailed spectra of SDSS J1216 with those 
of SDSSJ1035 in Fig. |SJ|, whi ch indicates t hat the orbital 
inclination is below about 75 ° jHellieill200ill . 

• The mass function of the orbit is (1.8 ±0.6) x 10"^ Mq. 
For a specified orbital inclination this gives the mass ratio 
of the stars. 

We have deliberately avoided using a mass-radius re- 
lation for properties of the secondary sta r. Empirical rela- 
tions exist for the mass donors in CVs (e.g. lSmith fc DhiUonl 
Il998h but the intrinsic complexity of these systems means 
that they are based on limited observational data. Theoreti- 
cal relations for CVs are not trustworthy at present because 
of the discrepancy between the predicted and observed prop- 
erties of the CV population (see Sect.0. Several empirical 
and theoretical relations exist for late-type main sequence 
stars which have evolved essentially independently, but re- 
cent detailed investigations of detached eclipsing binary sys- 
tems consistently find that the radii of late-type dwarfs are 
up to 10% greater than predicted iTorres &: Riba3 l2002l : 
iLopez-Morales fc Ri bas"'2005': lRibajl2005l) . 

The presence of double-peaked emission suggests that 
the o rbital inclinati on, i, of the SDSS J1216 system is quite 
high. lWarneil lll976l) noted that double-peaked emission lines 
indicate a CV with i > 50°, although there is no ob- 
vious reason why double-peaked li nes could not be see n 
from much lower inclination systems iHorne fc Marshligiil . 
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Choosing an inclination of 60° , using our value for the mass 
function, and assuming a mass of 0.6 Mq (1.0 Mq) for the 
white dwarf, we find that the secondary star has a mass 
of 0.022 ± 0.005 Mq (0.031 ± 0.005 Mq). These values are 
both in reasonable agreement with theoretical predictions 
of the properties of C Vs which have evolved past the or- 
bital period minimum jKolb fc Baraffelll999l : iHowell et all 
I2OOII) . Adopting a higher inclination would cause the in- 
ferred mass of the secondary star t o decrease. Finally , using 
MwD = 0.69 ± 0.13 Mq, as given bvl Smith fc DhillonI (Il998l) 
for CVs below the period gap, Fig. ll3l indicates that the sec- 
ondary star of SDSS J1216 has a mass in the range 0.012 to 
0.040 M0 (including the uncertainty in our measurement of 
the mass function). 

Fig. ll3l shows the effect of the above constraints on the 
masses of the components of SDSS J1216 and the inclination 
of the system. Light shaded areas indicate those parts of the 
diagram which are only accessible for an inclination of 75° 
or more. Dark shaded areas indicate the areas forbidden by 
the constraint on the maximum Keplerian velocity of the 
disc. It also shows the constraint due to the mass function 
for a number of selected inclinations between 20° and 75°. 

From Fig. I13l it is clear that there are two possible solu- 
tions to the system. The most likely of these, given that the 
double-peaked lines suggest a high inclination, is a donor 
star with a mass significantly below the hydrogen-burning 
limit. The less likely alternative is a white dwarf with a large 
mass of 1.15 Mq or greater. Thus SDSS J1216 is a good can- 
didate for a CV containing a brown-dwarf-mass donor star. 

SDSS J1216 resembles CD 552, which also has double- 
peaked lines, with a peak separation smaller than 
SDSS J12 16 but a very simi l ar ve locity amplitude of 17 ± 
4kms~ 'Hessman fc Hopd ll990l) . Based on their data, 
iHessman & Hopd ( 1990) suggested that CD 552 is a low- 
inclination CV containing a high-mass white dwarf and a 
main sequence secondary star. In the light of the predic- 
tions of theoretical models fe.g. lKolbHl993t iKolb fc Baraffd 
ll999l:lHowell et al.l200lt|Politanj2004ft . a viable alternative 
hypothesis is that CD 552 and SDSSJ1216 are CVs with 
substellar donor stars. This type of system is expected to 
be common and in fact should dominate the galactic CV 
population. We therefore feel confident in proposing that 
SDSS J1216 is a strong candidate for a CV with a brown- 
dwarf-mass secondary star. 



5 DISCUSSION AND CONCLUSIONS 

This work presents phase-resolved spectroscopy of six CVs 
recently discovered in SDSS spectroscopic survey observa- 
tions. For four of these systems we have obtained a pre- 
cise orbital period using radial velocities measured from the 
wings of H/3 or Ha emission lines. Three of the four systems 
have periods below the period gap: 96min for SDSS J0233, 
82 min for SDSS J1035 and 99 min for SDSS J1216 (although 
we have not been able to completely rule out the possibility 
that the 92 min one-day alias is in fact the correct period). 
The fourth, SDSSJ0911, has a period of 296 min, which is 
substantially longwards of the period gap. 

Whilst we have found a clear 96 min periodicity in the 
emission-line radial velocities of SDSS J0233, a period of 
60 min is present in the light curves of this system. We sug- 



gest that this system is a low-accretion-rate intermediate 
polar and the 60 min variation represents the spin period of 
the white dwarf. 

The spectra of SDSSJ1035 show it to be an eclipsing 
system; subsequent ULTRACAM photometry has found the 
eclipses to be total (Littlefair et al., submitted). Its short 
period and spectrum dominated by a cool white dwarf are 
reminiscent of WZ Sge systems (e.g. IPatterson et al.ll2005al : 
iTempleton et all 120061) . From modelling the spectrum we 
find a white dwarf effective temperature of 11 700 ± 1500 K 
and an upper limit of LO to the spectral type of the sec- 
ondary star, which is not detected in the spectrum. The 
high quality of the VLT/FORS2 observations has also al- 
lowed Doppler tomograms to be constructed for its accre- 
tion disc, and these show a clear bright spot and a disc with 
comparatively high-velocity motions. 

We have found SDSSJ1216 to have a very low veloc- 
ity amplitude of 13.8 ± 1.6kms~^. To our knowledge this 
is the smallest amplitude so far measured for a CV from 
emission-line radial velocities (excluding lower limits set by 
non-detection of velocity variation). Through applying a se- 
ries of physically reasonable constraints we find that the sys- 
tem must have either a high-mass white dwarf primary com- 
ponent and/or a secondary star with a mass in the brown 
dwarf regime. Systems with donor masses well below 0.1 M0 
are expected to be very faint, because the secondary star 
is dim and the mass transfer rate is low, and so under- 
represented in the current observational population of CVs. 

There are several candidates for CVs with brown- 
dwarf-mass secondary sta rs and their properties have 
been reviewed by iLittlefair et al.l ll2003ll . to whose list 
should be added the r ece ntly-s tudied CVs HS 2331-1-3905 
fArauio-Betancor et al.' l2005l) and RE J1255-I-266 
(Patterson et al. 2005b). In general, inferences of the 
secondary mass have been drawn from the spectral energy 
distribution of the system, from calibrations between super- 
hump period excess and mass ratio, or through detection of 
spectral lines of the s econdary component in a few cases. 
ILittlefair ot al.j i2003l) found that direct evidence for a 
brown dwarf secondary does not exist for any CV, but that 
some systems have substantial indirect evidence. 

Three of the four CVs studied here have very short or- 
bital periods, which is consistent with the expectation that 
studying successively fainter samples of CVs means that the 
targets become on average intrinsically fainter and of shorter 
period dPatterson et al.l2005bl:lGansick(i2005h . Our ongoing 
research program of CVs discovered by the SDSS will give 
orbital periods of a larger homogeneous sample of these sys- 
tems, allowing us to accurately assess the biases affecting 
the observed sample of CVs and so accurately measure the 
properties of the intrinsic population of these objects. 
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